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Abstract 

We study the coincidence problem of late cosmic acceleration by assuming that the present ratio 



energy component we consider two different candidates, first a quintessence scalar field, and then a 
tachyon field. In either cases analytical solutions for the scale factor, the field and the potential are 



between dark matter and dark energy is a slowly varying function of the scale factor. As dark 
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derived. Both models show a good fit to the recent magnitude-redshift supernovae data. However, 
the likelihood contours disfavor the tachyon field model as it seems to prefer a excessively high 

i ■ 

Q ■ value for the matter component. 
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I. INTRODUCTION 



Nowadays it is widely accepted that the present stage of cosmic expansion is accelerated 
albeit there are rather divergent proposals about the mechanism behind this acceler- 
ation. A cosmological model of present acceleration should not only fit the high redshift 
supernovae data, the cosmic microwave background anisotropy spectrum and safely pass 
other tests, it must solve the coincidence problem as well, namely "why the Universe is 
accelerating just now?" , or in the realm of Einstein gravity "why are the densities of matter 
and dark energy of precisely the same order today?" [3J] -note that these two energies scale 
differently with redshift. While it might happen that this coincidence is just a "coincidence" 
-and as such no explanation is to be found- we believe models that fail to account for this 
cannot be regarded as satisfactory. 

In a class of models designed to solve this problem the dark energy density "tracks" 
the matter energy density for most of the history of the Universe, and overcomes it only 
recently (see, e.g., Ref. j^J). However, these models suffer the drawback of fine-tuning 
the initial conditions whereupon they are not, after all, much better than the conventional 
"concordance" model which rests on a mixture of matter and a fine-tuned cosmological 
octant 0. 

There is an especially successful subset of models based on an interaction between dark 
energy and cold matter (i.e., dust) such that the ratio r of the corresponding energy densities 



tends to a constant of order unity at late times B, y, y, h h h thus solving the 

problem. However, the current observational information does not necessarily imply that r 
ought to be strictly constant today. For the coincidence problem to be addressed a softer 
condition may suffice, namely that at present r should be a slowly varying function of the 
scale factor with r(a = ao) ~ 3/7, the currently observed ratio. By slowly varying we mean 
that the current rate of variation of r(a) should be no much larger than H , where H = a/a 
denotes the Hubble factor of the Friedamnn-Lemaitre-Robertson- Walker (FLRW) metric, 
and a zero subscript means present time. It should be noted that because the nature of dark 
matter and dark energy is largely unknown an interaction between both cannot be excluded 
a priori. In fact, the possibility has been suggested from a variety of angles 

To avoid a possible conflict with observational constraints on long-range forces Q| we 
consider that the baryon component of the matter does not participate in the the interaction 
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and, further, to simplify the analysis -i.e., in order not to have an uncoupled component- 
we exclude the baryons altogether. While this might be seen as a radical step it should 
be taken into account that our study restricts itself to times near the present time and 
these are characterized, among other things, by a low value of the baryon energy density 
(5% or less of the total energy budget, approximately six times lower than the dark matter 
contribution and fourteen below the dark energy component whereby it should not 
significantly affect our results. This is in keeping with the findings of Majerotto et al. [16]. 
For interacting models encompassing most the Universe history in which the baryons enter 
the dynamical equations as a non-interacting component, see Refs. Q and Q. 

The target of this paper is to present two models of late acceleration that fulfill "soft 
coincidence" , namely: (i) when the dark energy is a quintessence scalar field, and (ii) when 
the dark energy is a tachyon field. The latter was introduced by Sen Q| and soon afterwards 
it became a candidate for driving inflation as well as late acceleration -see e.g., Ref. jit] ]. 
The outline of this paper is as follows: Section II considers the quintessence model with a 
constant equation of state parameter. There it is assumed that the quintessence field slowly 
decays into dark matter with the equation of state of dust. Section III considers the tachyon 
field and again assumes a slowly decay into dust. This time, however, the equation of state 
parameter is allowed to vary. Finally, section IV summarizes our findings. 

II. THE QUINTESSENCE INTERACTING MODEL 

We consider a two-component system, namely, cold dark matter, described by an energy 
density p m , and a quintessence scalar field with energy density and pressure defined by 

P0 = |0 2 + W), and p<p = l^-V^), (1) 

respectively, in a spatially flat FLRW universe. The over dot indicates derivative with 
respect to the cosmic time and V(<f>) is the quintessence scalar potential. We assume that 
these two components do not evolve separately but interact through a source (loss) term 
(say, Q) that enters the energy balances 

p m + 3Hp m = Q (2) 
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and 



p^ + 3H(p^, + pf) = -Q, 



(3) 



where in view of (JTj) last equation is equivalent to 



+ 3H<j) + V 



-Q. 



(4) 



In the following we constrain the interaction Q by demanding that the solution to Eqs. 
(j2J) and (JHJ) be compatible with a variable ratio between the energy densities r(x) = p m / P</>, 
where x = a/ao is the normalized scale factor, and that around the present time r(x) is a 
smooth, nearly constant function with r(x = 1) = tq being of order one. We also assume 
that the quintessence component obeys a barotropic equation of state, it is to say = w^p^ 
with a negative constant (a distinguishing feature of dark energy fields -quintessence fields 
or whatever- is a high negative pressure). In virtue of these relations the set of dynamical 
equations reduces to a single equation 

P0 + 3# (l + ^^)p, = 0, (5) 



whose solution is 

,(0) -3I(x) 



p+{x)=pfe*<™, (6) 



with 

f x , ,x dx' , „ Wj, + hxr'(x) , . 

I( x )= F(x')— , where F(x) = 1 + ^ 3 1 J , 7 
Ji x' 1 + r(x) 

and the prime means derivation with respect to x. On the other hand, by combining Fried- 
mann's equation 
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3H 2 = K(p m + p$) (k = 8nG) 



(8) 



with Eq. (0) we get 



m*) = H.ii±g&6-t™, ( 9) 



where Hq = y p^K,/3 denotes the current value of the Hubble factor. From this, it follows 
that 

If this integral could be solved analytically, we would obtain the scale factor in terms of the 
cosmological time. 

Equations (J3J) and © alongside (JBJ), ©, and © imply 



Q(x) - 



(o) 



V 1 + r o \/l + r(x) 



xr'ix) , , 

— w^r{x) 



e-pW. (11) 



Note that Q(x) is a positive-semidefinite function, as it should. A negative Q(x) would imply 
a transfer of energy from the matter to the scalar field which migth violate the second law 
of thermodynamics. While in view of the unknown nature of dark matter and dark energy 
we cannot say for certain that these components fulfill the aforesaid law, in the absence of 
any evidence against it, the most natural thing is to assume that they obey it. 

From the definitions ((H) and the equation of state = w^p^ the quintessence field and 
its potential are given by 



/3(l + iiu) f x 1 dx' 



and 



V(x) = V e- 3I{x) , (13) 
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respectively. Here 0o is an integration constant, and Vq = \ {1 — w 4 



As said before, we apply the above formalism to the case in which the variable x is not 
far from unity whence the ratio r(x) can be approximated by 

r(x) ~r + e (l-x), (14) 



where r is the present value of the ratio between the energy densities p m and p^, and Eq is 
a small positive-definite constant. We do not consider negative values for Eq since it would 
imply that r(x) was increasing in the recent past and therefore that it oscillates. While 
we are unaware of any definitive argument against this possibility it looks certain that only 
contrived models may lead to this behavior. Further, oscillations in r(x) may seriously 
jeopardize the well tested picture of structure formation j^j]. (Note in passing that the 
choice ffTlJl implies that | r(x) |< H for x ~ 0(1)). 



It follows that 



7-1/ \ , , Oil ~ & 2 X nK s 

F(x) = l-\ , (15) 

a 3 — a 4 x 



where the a, are constants given by 

a 1 = w (/) , a 2 = ho, a 3 = l + r + £ , a 4 = e , (16) 



respectively. With this, we obtain 



Pm{x)=p\'\ } (a 3 -l - otiXjx \ <*» /. (17) 

«3 — ci±x 



and 



«3 — OI4X 



P4>{ X )=PV[ — Z ) x V 3 ; > ( 18 ) 
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The Hubble parameter can be written as 



3 / °l+°3 \ 

H(x) = H x n - 3 J, (19) 



where Hq was given above. From this expression the scale factor is shown to follow a power- 
law dependence on time 



, t \ 3 ( Q l+ a 3) 

a(t) = «o - • (20) 



For the Universe to accelerate, the constraint „, 2 " 3 ? > 1 must be fulfilled, i.e., 

' 3(01+03) ' ' 

l + r + £ <-3^. (21) 

This, alongside the condition that the energy densities decrease with expansion implies 

-W4 < 1 + r + e < -3itty. (22) 
Fig. [T] shows a good fit to the "gold" set of supernovae data points of Riess et al. Q|. 

The likelihood contours are depicted in Fig. El The mean values of the free parameters 
are: = 0.78538, = -0.757284, e = 0.0777764. Notice that fi is above the con- 
cordance ACDM value of 0.73 and that is significantly larger than the value found in 
non-interacting models. However, we have not considered phantom fields (scalar fields as 
given by Eq. (0) do not encompass phantom behavior), otherwise a shift of toward more 
negative values should be expected. Notice (top panels) that the parameter Eq is rather 
degenerate. 

Under restriction (J2TJ) the interaction term reads 



Q(x) = -3H p ( ; 



(o) 



3[(«3 — 1 — a.±x)\ 



(03-1-042;) -§ f 21 ^) , OQ> 
(a 3 -a 4 )— -2-ac V 3 /, (23) 

(«3 — 04a;) 
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FIG. 1: Distance moduli vs redshift for the quintessence-dark matter interacting model. In plotting 
the graphs the expression \x = 51og(ii + 25, with di = (1 + z) Jq H~ l (z')dz' in units of megaparsecs 
was assumed. Here = —0.75, eo = 0.077 and tq = 0.27. For comparison we have also plotted 
the prediction of the concordance ACDM model with ft, m o = 0.3. The data points correspond to 
the "gold" sample of type la supernovae of Ref . |3] • 

see Fig. El 

In its turn, the scalar field <fr and the scalar potential V(x) are given by 



(x) = O 



tanh 



-i 



tanh 



-i 



«3 



L x 



I — 2A 

as 



(24) 



and 



V(x) = V I x V • , 

a 3 — a^x 



(25) 
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FIG. 2: Likelihood contours for the quintessence-matter interacting model showing the 68%, 90% 
and 99% confidence intervals. The likelihoods are marginalized over the rest of parameters. The 
prior £l m + = 1 was used. The right bottom panel shows the probability function of the 
quintessence parameter density. 



respectively. Here </> 



12 / 1+Wd 



j tanh 1 


Jl - ^ 




.V Q3 . 



and Vo = |(1 - w^)p^\ see Fig. H 



Equations ((23) and ((23) lead to 
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Figures E] and El taken together show that the potential decreases with the Universe expan- 
sion. While we do not know whether there is any field theory backing this potential it is 
intriguing to see that around = it behaves as 



V(</>) ^C l( j)- 2 + C 2 + C 3 <f) 2 + C 4( f) A + ... (27) 



where the CVs are constants. The first term is used in quitessence models -see e.g., Ref. |21[, 
whereas the third and fourth terms of the expansion are well-known potentials in inflation 
theory (chaotic potentials) |22j; the second term plays the role of a cosmological constant. 
This leads us to surmise that, in reality, V(<fr) might be considered an effective potential 
resulting from the combination of a number of fields. 

III. THE TACHYON INTERACTING MODEL 

The tachyon field naturally emerges as a straightforward generalization of the Lagrangian 
of the relativistic particle much in the same way as the scalar field <fi arises from generalizing 



L " sr " ,is ' 1 r " 1 "■• 



tiality as dark matter |24j and dark energy |25( has awakened the interest in it. We begin 
by succinctly recalling the basic equations of the tachyon field to be used below where its 
interaction with cold dark matter (dust) will be considered. 
The stress-energy tensor of the tachyon field 

= / . V } (P \ [-9a b (1 + <^V,c) + V,a<P, b ] , (28) 
admits to be written in the perfect fluid form 

T ab } = Pv U °. U b + P<p {9ab + U a U b ) , (29) 
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where the energy density and pressure are given by 



V 1 - <P 

respectively, with 



n , = , and Pv = -V{<f>) y/1 - (30) 



ip = ip a u a = J -g ab (p a (p h , and u a = —r L , where u a u a = -l. (31) 
v ip 

In the absence of interactions other than gravity the evolution of the energy density is 
governed by p^ = —3H(p 2 p<p, therefore when ip 2 < 1 it decays at a lower rate than that for 
dust. It approaches the behavior of dust for ip 2 — > 1, thereby in this limit the tachyon field 
behaves dynamically as pressureless matter does. Consequently, we shall assume up 2 < 1 
since for <p 2 = 1 both components obey the same equation of state for dust. 

For an interacting mixture of a tachyon field and cold dark matter, with energy density p m 
and negligible pressure, the interaction term, Q, between these two components is described 
by the following balance equations 

p m + 3H(p m + II m ) = Q , (32) 



Pip + 3E(p 2 p^ = -Q . (33) 

The Il m term, on the left hand side of Eq. (JH2J), accounts for the fact that the matter 
component may be endowed with a viscous pressure or perhaps it is slowly decaying into 

In either case one can model this term as Il m = ap m H 



dark matter and/or radiation 
with a a small negative constant since Il m is a small correction to the matter pressure -see 
[22! and references therein. 

As before, we consider the ratio between the densities of matter and tachyonic energy a 
function r(x) of the normalized scale factor (to be specified later), and again, we must have 
Q(x) > -its expression is to be found below. Then, equations (J3~2"j) and (fBlIj) combine to 
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The latter can be solved to 



P<p(x) = P^e~ 3Hx) , where I(x) = / F(x') % 

J i x 



(35) 



with 



F(x) = 1 + 



w,p(x) + ar(x) + ^r'(x)x 
1 + r(x) 



The interaction term takes the form 



(36) 



Q(x) = 3p<pH[ 



x) 



r(x) 



1 + r(x) 



a - w^x) + 



xr'(x) 
3r(x) 



-31 (x) 



(37) 



and the tachyonic scalar field and its potential obey 



ip(x) = (p + 



VI + r f x j l + w {p (x') s/^/) dx' 



H J 1 V l+r{x') 



x 



(38) 



and 



V(x) = V J-w v (x) e~ 3/ >\ 



(39) 



respectively. 

Up to now we have left the ratio function r(x) free. As before, we specify it for x values 
around unity as 



r(x) ~ r + £q(1 - x ): 



(40) 



where r = (p m / P<p)o, an d £ o is once again a small positive-definite constant. Likewise, we 
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assume that the equation of state parameter w v is given by 



w v (x) =w + wi(1 - x), 



(41) 



where wo and w± are constants, the first one denotes the current value of the w^x) function, 
and the second one is minus its first derivative, which is expected to be small. Thus, 



F(x) = 1 + 



ai — b\X 
a 2 - b 2 x ' 



(42) 



where the constants a,- and b, } stand for 



ai = ar + wq + wi + e , a 2 = 1 + r + e 



(43) 



and 



bi = wi + (a + l)e , b 2 = e , 



(44) 



respectively. It follows that 



p m {x) = pf (1 + r ) 3ft x' 3132 [r + e (l - x)] [1 + r + e (l - x)Y 



3/3i 



(45) 



as well as 



p v (x) = pf (1 + r ) 3ft a;- 3 ^ 2 [1 + r + e (l - x)]~^\ 



(46) 



with 



and 



°2 a 2 



l + r + e 



ai l + wi + 2e + r (l + a) 

P2 — J- H — : 

a2 1 + r + e 



The Hubble function 



3/3! -1 -302 (-3/3i + l) 

if (a) = ff (1 + r )^^ x— [1 + r + e (l - x)] 3 , 



(47) 
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follows from the Friedmann's equation. 

Although last expression is comparatively simple, the scale factor derived from it is not 



^(1 + ro)^ H (t - t ) = + r + e (l - x))^ 



x 2 F X 



3& 1-3A 
2 ' 2 



n 



e x 



1 + r + e 



Ci 



where 2 -Pi is the hypergeometric function 28] and 



(48) 



Ci 



1 + ro 1 



e 



1 + r + e 



2 F 1 (x = l). 



Fig. |H1 portrays the evolution of the scale factor in terms of the cosmological time as well 
as the deceleration factor q = —a/(aH 2 ) versus the redshift for two selected values of the 
parameters. 



As Fig. [7| shows the model fits the supernova data points not less well than the con- 
cordance ACDM model does. The likelihood contours, Figs. [HI and El were calculated with 
the method of Markov's chains. We used the prior Q m + Q v = 1 and that the parameters 
Wq and w\ are restricted by the condition that the value of the right hand side of Eq. (14 1|) 
must lay in the interval [—1, —1/3). The mean values of the parameters are: Vt^ = 0.246, 
wo = —0.773, wi = 0.22, e = 0.0087, a = —0.76. Here, e$ is not so weakly constrained 
by the supernovae data as in the quintessence model. The present model predicts a mild 
evolution of the equation of state parameter with redshift. This is slightly at variance with 
the findings of Jassal et al. |30j , but agrees with the model independent analysis of Alam et 
al. bJ. 



The interaction term is given by 

(1 + ro)^ 

Q(x) = Qo — 7 r 7 {3 a -w - wx(l - x) r + e (l - x)\ - e x} 

[3r (a - w ) - e \ 
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-9/32 

XX 2 



{1 + r + e (l -x)} 



1-901 



(49) 



with Q = \p\> H [3r (a - w ) - e ]. 

Likewise, the tachyon field and the potential are found to be 




o Jl 



and 



V(x) = pf(l + r ) 3A yZ-wo-w^l-x) x" 3/32 [1 + ro + e (l - x)]- 



(51) 



respectively -see Fig. [TUl 

IV. CONCLUDING REMARKS 

We have studied two models of late acceleration by assuming (i) that dark energy and 
non-relativistic dark matter do not conserve separately but the former decays into the latter 
as the Universe expands, and (ii) that the present ratio of the dark matter density to dark 
energy density varies slowly with time, i.e., | r |o< Ho. This second assumption is key to 
determine the interaction Q between both components. 

In the quintessence model (section II) we have considered the equation of state parameter 
constant while in the tachyon field model (section III) we have allowed it to vary slightly. 
Actually, there is no compelling reason to impose that this parameter should be a constant. 
However, Jassal et al. |30j have pointed out that the WMAP data [15| imply that in any 
case it cannot vary much. By contrast, Alam et al. using the sample of "gold" supernovae 
of Riess et al. find a clear evolution of w in the redshift interval < z < 1; however when 
strong priors on Q m0 and H Q are imposed this result weakens. Nevertheless, the analysis of 
these two papers assume that the two main components (matter and dark energy) do not 
interact with each other except gravitationally. The parameter u> presents degeneration in 
both models, therefore we must wait for further and more accurate SNIa data, perhaps from 
the future SNAP satellite, or to resort to complementary observations of the CMB. 

In both cases (quintessence and tachyon), we have found analytical expressions for the 
relevant quantities (i.e., the scale factor, the field and the potential) and the solutions 
are seen to successfully pass the magnitude-redshift supernovae test -see Figs. [T] and [7| 
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Nevertheless, it is apparent that the the tachyon model favors rather high values of the 
matter density parameter (see bottom right panel of Fig. |HJ) which is at variance with a 
variety of measurements of matter abundance at cosmic scales 32j which, taken as a whole, 
hint that Q m should not exceed ~ 0.45. In consequence, the quintessence model appears 
favored over the tachyon model. Our work may serve to build more sophisticated models 
aimed to simultaneously account for the present acceleration and the coincidence problem. 

Previous studies of interacting dark energy aimed to solve the coincidence problem by 
demanding that the ratio r be strictly constant at late times needed to prove the stability 
of r at such times. This was achieved by showing that the models satisfied an attractor 
condition that involved the equation of state parameter of matter and dark energy as well 
as the Hubble factor and its temporal derivative Since in the case at hand the 

coincidence problem is solved with a (slowly) varying ratio r no stability proof is necessary 
at all and no attractor condition is needed. 

Our analysis was confined to times not far from the present (i.e., for x ~ 0(1). To recover 
the evolution of the Universe at earlier times (x <C 1), when the matter density dominated 
and produced via gravitational instability the cosmic structures we observe today, we must 
generalize our study along the lines of Refs. B and Q and include the baryon component 
in the dynamic equations as an uncoupled fluid. 

We restricted ourselves to scenarios satisfying w > —1. Scenarios with w < — 1 (the 
so-called "phantom" energy models) violate the dominant ener gy c ondition though, never- 
theless, they are observationally favored rather than excluded |3J| and exhibit interesting 
features |24| that might call for "new physics". We defer the study of phantom models 
presenting soft coincidence to a future publication. 
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FIG. 3: Evolution of the interaction term with expansion. The lower eq, the slower the decrease 
of the quintessence field. Here = —0.8 and ro = 3/7 
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FIG. 4: Evolution of the scalar field with expansion. Again, = —0.8 and ro = 3/7. 
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FIG. 5: The effective potential versus the scalar field. Here w<$, = —0.8 and ro = 3/7. 
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FIG. 6: The left panel shows the scalar factor as a function of cosmological time. The right panel 
shows the deceleration parameter q as a function of the redshift z = x~ l — 1. In both cases u>o=-0.9, 
wi=0.002, r = 3/7 and a=-0.2. 
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FIG. 7: Distance moduli vs redshift for the tachyon-dark matter interacting model. In plotting the 
graphs the expression \x = 51og(iL + 25, with di, = (1 + z) J Q H~ 1 (z')dz' in units of megaparsecs 
was assumed. We have taken the values of the best fit model, namely: wq = —0.99, w\ = 0.95, 
£o = 0.0042, a = —0.98 and ro = 0.136. For comparison we have also plotted the prediction of the 
concordance ACDM model with f2 m o = 0.3. The data points correspond to the "gold" sample of 
type la supernovae of Ref . [3] • 
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FIG. 8: Likelihood contours for the tachyon-matter interacting model (w x ,o vs. eo -top panel-, and 
Q.tp vs. w x ,o, bottom panel) showing the 68% and 98% confidence intervals. The likelihoods are 
marginalized over the rest of parameters. 
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FIG. 9: Top panel: likelihood contours for the tachyon-matter interacting model (a vs. £l m )) 
showing the 68% and 98% confidence intervals. Bottom panel: probability function of the matter 
density parameter. The likelihoods are marginalized over the rest of parameters. 
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FIG. 10: The effective potential as a function of the tachyonic field. Again, wq=-0.9, w\= 0.002, 
r = 3/7 and a= - 0.2. 
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